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Abstract. The scalar-tensor Dirac-Born-Infeld (DBI) inflation scenario provides a simple 
mechanism to reduce the large values of the boost factor associated with single field models 
with DBI action, whilst still being able to drive 60 efolds of inflation. Using a slow-roll 
approach, we obtain an analytical expression for the spectral index of the perturbations and, 
moreover, determine numerically the regions of the parameter space of the model capable 
of giving rise to a power spectrum with amplitude and spectral index within the observed 
bounds. We find that regions that exhibit significant DBI effects throughout the inflationary 
period can be discarded by virtue of a blue-tilted spectral index, however, there are a number 
of viable cases — associated with a more red-tilted spectral index — for which the boost 
factor is initially suppressed by the effect of the coupling between the fields, but increases 
later to moderate values. 
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1 Introduction 

Models of inflation based upon the Dirac-Born-Infeld (DBI) action [1, 2] have been much 
studied in the past few years as an well-motivated example of string cosmology. In this 
context, the field x appearing in the DBI Lagrangian is directly related to the radial coor- 
dinate of a D3 brane moving in a 'throat' region of a compactified space with a speed limit 
imposed upon its motion. The DBI model is an interesting example of k-inflation [3] (in 
which the Lagrangian is an arbitrary function of the inflaton and its kinetic term) and has 
been much studied in this regard. K-inflationary models are noteworthy as the sound speed 
of the perturbations Cg can be less than 1, which has the effect of forcing the scalar per- 
turbations generated during inflation to freeze-in at scales shorter than the curvature radius 
[4]. In the DBI model the sound speed is given in terms of the boost factor 7 (which plays 
a role analogous to the Lorentz factor in special relatively) by = 7"^, so c<j — )• as the 
speed limit is saturated. The original DBI papers concentrated mainly on this 'relativistic' 
regime where 7 is large and for which a simple power law inflation solution can be obtained. 
However, large values of the boost factor are inextricably linked with excessive amounts of 
non-Gaussianity in the bispectrum of the CMB fluctuations as the non-linearity parameter 
/nl is proportional to 7^ in the relativistic case. This leads to an upper bound on the value 
of the boost factor 7 < 20, assumed to apply during the first 10 of the last 60 efolds of 
inflation which directly affect the CMB radiation. Since the boost factor increases with the 
number of efolds of inflation A^maxj many DBI models with A'max sufficiently large predict 
a high level of non-Gaussianity exceeding current observational bounds (cf. [5, 6] and also 
[2, 7]). As shown in [8] (cf. [9]) including non-Gaussianity constraints gives an upper bound 
on the slow-roll parameter. This greatly restricts the number of efolds of inflation driven by 
the DBI field. 

One should also take into account the theoretical bound on the range of Xi which arises 
simply because the throat must have finite size (sometimes referred to as the bulk-volume 
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bound). This is related to the level of tensor perturbations (given in terms of the the tensor- 
to-scalar ratio r) from inflation via the Lyth bound, 



which is true for both DBI and slow-roll inflation [7, 10]. This gives an upper bound on r 
that can be combined with the lower bound derived from the upper bound on 7 to place 
severe constraints on ultra-relativistic DBI models. Similar ideas are used in [11] to derive 
a series of consistency relations between the principle observable quantities. The problems 
with observational constraints are made explicit in two detailed numerical studies [5] and [6] , 
which use Monte Carlo methods to compare single field DBI models to recent cosmological 
observations. Both studies find that the majority of models cannot satisfy the bulk- volume 
bound; those that remain are slow-roll DBI models with small values of 7. 

Several routes have been explored that have a bearing on the problem, in particular the 
inclusion of other fields in the determinant in the DBI Lagrangian to represent non-radial 
directions in the throat. With extra degrees of freedom, however, come extra features, in par- 
ticular non-radial trajectories of the probe brane and non-adiabatic (entropic) fluctuations. 
These were addressed in [12] (see [13] for recent developments), where the authors introduced 
'spinflation', accounting for angular momentum in the UV DBI model. Other models involv- 
ing the angular coordinates of the DBI scenario have been considered in [14, 15] (see also 
[16-20]). A different type of multi-field model was introduced in [21] in which inflation is 
driven by two standard DBI fields, each corresponding to a brane with its own sound speed 
(see also [22], which generalises the assisted inflation scenario for multi- field DBI inflation). 
Another way in which the standard DBI scenario can be modified is by coupling the DBI 
field to another scalar field, or, as in [23, 24], to the gravitational term in the action. The 
effect of coupling in multifield DBI models was considered in [25, 26], focusing on the effect 
of particle production due the interaction between the 'inflaton' brane and trapped branes 
in the warped throat.^ 

In a previous work [27], we studied the consequences of coupling the DBI action to 
a canonical scalar field in a scalar-tensor theory, in order to investigate the effect of an 
additional scalar field with a non-minimal coupling to gravity in the effective action. We 
studied a coupling of the form A = exp{(3(p), arising due to a conformal transformation to 
the Einstein frame. As in the two-field model described in [28], when the coupling is non-zero, 
one field acquires a large effective mass, even for couplings of order a tenth of the strength of 
gravity (/? > 0.1). The parameters of the DBI field, which are dependent on the additional 
scalar field, vary during inflation so that the number of e-folds is extended and the boost 
factor is decreased (when compared to standard DBI inflation with the same bare param- 
eters). Two potentials for the canonical field were investigated, one with a minimum and 
one without, yielding similar results for both the background and perturbations, suggesting 
that the conclusions hold for any choice of potential steep enough to allow ip to closely track 
the minimum of its effective potential. As discussed above, many DBI models predict a high 

^ Although, hke the scalar-tensor DBI inflation model in [27], this scenario involves a DBI action non- 
minimally coupled to a second field, the motivation and dynamics are very difi'erent: for example, in our 
model the energy density of the yi-field is not negligible (and indeed, is often the dominant contribution to 
the right-hand side of the Friedmann equation) and varies slowly as the DBI field evolves. These differ- 
ences notwithstanding, something like this scenario might perhaps provide an promising framework for the 
phenomenological model discussed here. 
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level of non-Gaussianity arising due to large values of 7: in the coupled model, the boost 
factor is reduced, so the level of non-Gaussianity (which is proportional to 7^ in single field 
DBI models) may be expected to be smaller than the standard DBI case. If this is the case, 
the presence of the canonical scalar field may alleviate some of the problems of the DBI 
inflationary scenario, although we do not address the issue of non-Gaussianity in this paper. 

To understand the significance of the model, it is essential to calculate the predicted 
values of any observable parameters so as to compare with data, however, this is made rather 
difficult in this case as there are a large number of free parameters capable of affecting the 
results. In this article, we approach the problem from two fronts. After a brief resume of the 
particulars of the scalar-tensor DBI model in Sec. 2, we focus on the 'slow-roll' limit of the 
model in Sec. 3, showing that with a judicious choice of slow-roll parameters, the equations 
of motion for the perturbations can be made tractable and solved (to first-order) to obtain 
an expression for the spectral index. In Sec. 4, it is shown that the parameter space can 
be severely constrained by considering the values of the power spectrum amplitude and the 
spectral index. The latter measurement in particular rules out a large class of models — 
those with small values of Cg at horizon crossing, which give rise to a blue-tilted spectrum 
— so that only scenarios in which DBI characteristics are suppressed at the beginning of the 
observable period of inflation are allowed. Finally, our results are summarised in Sec. 5 



2 Scalar-tensor DBI inflation 

In the Einstein frame (cf. [27] for details) the scalar field (p has a canonical kinetic term and 
does not couple explicitly to the gravity sector, so the action for the model takes the form 



S = j d^x^ [^R + C^ + £dbi] , (2.1) 



where 



£^ = --gt^-^^^^^^ - U{p), (2.2) 
£dbi = A\ip) [f-\x) (1 - 7-') - Vix)] . (2.3) 

^(x) = ^'^^X^ s-i^d U{ip) are the potentials for the DBI field x ^iid the canonical scalar field 
ip respectively. Two forms for U{(p) were investigated in [27], one of an exponential form 

U{p) = Uoexp{-rip), (2.4) 

and another, quadratic potential U{if) = Uo{(p — rj)'^ with a minimum at f = rj. It was found 
that when the coupling between the fields was non-zero, both potentials led to extremely 
similar behaviour, so for simplicity the exponential potential (2.4) will be used throughout 
the body of the article, and the results of the numerical analysis for the quadratic potential 
will be summarised in Appendix B. 

/(x) is the warp factor of the DBI field, a measure of the geometry of the throat region 
of the compactified space in which the brane moves. As in [27], in the following the 'mass-gap 
solution' [2, 29] 

will be used, which exhibits the salient features of a capped throat with a cutoff at the IR 
end. The presence of the mass scale (which must be subplanckian in order to get a warped 
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throat) suppresses the DBI effects as the X'^^ld moves to smaller values. Generally, the 
smaller the value of the larger the 7 factor can become during the inflationary period and 
the more efolds one observes. 

The boost factor 7 appearing in (2.3) is modified by a (/9-dependent factor of the con- 
formal coupling A{ip) compared to the standard DBI model by a factor, so that 

7= , ^ , (2.6) 

where the dot indicates a derivative with respect to cosmic time. (Here, as throughout this 
article, aflat Friedmann-Robertson- Walker metric is assumed: ds^ = —dt^ + a?'{t)5ijdx'^dx^.) 
The conformal coupling A{ip) is given by 

A{^) = exp(/3(^), (2.7) 

in terms of the constant /3. 

The dynamical evolution of the system is determined by the background equations of 
motion for the homogeneous fields 



X + 3Hj-'x + M ^1(1 - 37"' + 27"') + A'j-'V^ = -(3x^{3j-' - 1), (2.8) 

^ + 3Hip + U^ = /3A^ [f-\i - 37"^ - 7) - , (2.9) 
where fy^ = df /dx and = dV/dx, together with the Friedmann equations 

3H^ = l^^ + U + A^[r\-f-l) + V], (2.10) 

-2H = ^^+jA^X^. (2.11) 



3 Analytical treatment of the perturbations 

In [27], the equations of motion for the first-order perturbations and the perturbed Einstein 
equations were derived in full and solved numerically. As was expected, the perturbations of 
the DBI field were found to propagate with a sound speed defined by = 7"^; however, many 
extra terms appear involving factors of Cg, the coupling /3 and derivatives of the warp factor 
/, which render the full equations intractable. In this section, we derive approximations for 
the background dynamics in the scenario in which the (^-field is trapped in the minimum 
of its effective potential. It is shown that in this case one can rewrite the perturbation 
equations in terms of slow-roll parameters and hence obtain an expression for the spectral 
index of scalar perturbations that is in excellent agreement with the numerical analysis in 
the following section. 

3.1 Defining slow-roll parameters 

As was shown in [27], in the coupled case with /3 > 0, the situation is simplified considerably 
as the extra mass term on the right-hand side of (2.9) means that the effective potential can 
have a minimum value at which the field value <p = fmin satisfies, 

- f3e^l^'P^^^T^^j = 0, (3.1) 
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where T^^j = /~^(4 — 7 — 37~^) — AV « —4V is the trace of the 'bare' DBI stress-energy 
tensor and the ip dependence of A{ip) has been written exphcitly using (2.7). With the 
exponential potential (2.4) the condition for the minimum in this case is 

U ~ ^A^V, (3.2) 
V 

so the minimum is given by the logarithmic function 

(cf. [28], in which a similar condition is obtained for two coupled scalar fields). Even for 
relatively small values of the coupling /3, the (/3-field is forced into this minimum, the position 
of which changes as the x-field evolves as^ 

so that the second Friedmann equation can be approximated as 

- 2ij ~ ASx^. (3.5) 

When the coupling between the two fields is switched on, the DBI equation of motion is 
dominated by the Hubble friction and potential terms (we have checked this numerically and 
found it to be an excellent approximation). Thus (2.8) can be approximated by the slow-roll 
equation 

3H^-^X + ^'7"'^x - 0' (3-6) 

i.e. 

A'^V 

Finally, as the contribution of the kinetic terms of the fields in (2.10) is negligible, one can 
rewrite the Friedmann equation as 

3H^ ^{l + f) A^V, (3.8) 

using (3.2). We are now in a position to define the first slow roll parameter e = —H/H^, 
which using (3.7), can be written 

1 4«\"2/K^2 



2^2^ V ; \v 

We can use this expression to obtain an expression for ip in terms of e. Substituting (3.7) 
into (3.4) gives 

P ^ iHe. (3.10) 



This corrects a minor error in eqn. (3.31) of [27]. 



- 5 - 



The rate of change of the slow-roh parameters should be second order in e. Differentiating 
(3.9) with respect to time gives 



1 



2He 



2He 
2Hc, 



2/3 



1 



2Hcs A^-i 



1 + M 



1 + M 



-1 V, 



XX 

V 




1 , 



<5 + 2 1 + 



3/3 



where the slow-roll parameters 5 and s are defined as follows 



5 = 



1 



4/3 



■ 1 + - 



-1 K 



XX 

V ' 



He, 



(3.11) 



(3.12) 



5 is the equivalent of the commonly used r/ parameter and s is a measure of the rate of change 
of the sound speed, which matches the definition used in the DBI literature (cf. [7]). Note, 
that s has an implicit dependence of A{ip) through the definition of 7 in (2.6). 

3.2 Rewriting the equations of motion 

It was demonstrated in [27] that when the (/?-field is trapped in the minimum of its effective 
potential, the effect of the perturbations 5ip{t^ x) is negligible and the coupled DBI model 
is effectively a single field system. The first-order perturbations can then be well-described 
by the Fourier components of the single variable u = ry.Qy.. Here is the combination of 
background quantities 

ry. = aAc~'^l^ (3.13) 

and = 5x + ■^'^ is the gauge- invariant Mukhanov-Sasaki variable constructed from the 
DBI field perturbation 5x and the metric perturbation ^. The perturbation equation in 
conformal time r (with ' = d/dr) is therefore [cf. eqn. (4.22) in [27]] 



4' + 



^x. 



Uk = 0, 



(3.14) 



with 



a 



A'^Xfx 



XX 



(1 



H /2' 
+lA\l-Csf 



f 



A^ 
He, 



-X 



+ |ylVc7'(l + c^ 



V.4 



-A^eZ 



^Vil + c?)4K + ^(1 + cl) + cU-'Vr,,,. 



2m 



Am 



H 



(3.15) 



Analysing each term in (3.15) separately, one can rewrite the expression in terms of the 
slow-roll parameters e, 5 and s as 



"-"xx 



3F^ (5 - s - 2e 



1 + ^(1 



(3.16) 
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The interested reader is directed to Appendix A for full details of the derivation of this 
equation. Differentiating with respect to conformal time gives 



nii + 



n 



Cs-H 



nil + fe 



(3.17) 



Differentiating again (again assuming the derivatives of e and s are second order) yields 

2 



(1 
2 - 



1+2^6 



1 _ M 



3s 



3s 



(3.18) 



using %' = T-L^il — e). Combining (3.16) and (3.18) means the combination that appears in 
the equation of motion can be simplified to 



XX 



3<5+(5 + f (2 



Defining the background variable z by 



CsH 



CsH 



a7V2e, 



one can show [using (3.11) and (3.19)] that 

r" 



'XX 



(3.19) 



(3.20) 



(3.21) 



which is correct to first order in e, s and 6. The first term on the RHS appears (without 
the factor of A in the definition of z) in single field k-infiation models [4]. Alone, this would 
indicate that the combination \uk/z\ becomes approximately constant for small k] the small 
correction proportional to (1+c^) arises due to the fact that (3.14) is an approximation to the 
full, two-field model given by eqns (4.21) and (4.22) in [27]. The correction is proportional to 
e and thus subdominant to z" / z in the mass term: this can be seen when comparing the size 
of the terms in the full equations (when evaluated numerically, as in the following section) 
and in many cases the correction is not too much larger than the terms proportional to the 
scalar field perturbations that have been neglected in (3.14). Thus, to a good approximation, 
we can neglect the second term in (3.21). In this case the evolution of is determined by 
the (approximate) equation 

z 



Uk ^ 0, 



(3.22) 



or, in terms of the slow-roll parameters 



V J 



Uk ~ 0. 



(3.23) 
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The comoving curvature perturbation of the full two-field system is given by eqn. (4.27) in 
[27] as 

where and are the gauge-invariant Mukhanov-Sasaki variables for the perturbations 
of the X and if fields respectively. Neglecting the Qip term and using (3.20), one gets 

as = Uk/r^ = \f2ecj'^ A~^U}^I z. The power spectrum is given by 

2 



3.3 Solving the perturbation equation 



z 



(3.24) 



As in the standard slow-roll case, eqn. (3.23) admits a solution in terms of Hankel functions. 
However, to account for the fact that Cs is not constant we follow the approach of [30] and 
rewrite the equation using the time variable 



y = 



Id! ~ 



(3.25) 



so ?/ = 1 at sound horizon crossing. For a related approach, see [9, 31]. The derivatives of 
can be rewritten in terms of the slow-roll parameters as 



duk 
dT 



-Csk{l — e — s) 



duk 
dy 



(3.26) 



(using n' = n'^il - e)) and 



d Uk _ 2 



[l-e-s)V 



d Uk 
dy"^ 



s(l - e - s)y 



duk 
dy 



(3.27) 



Substituting these expressions into (3.19) gives 



/, \2 2^^^fc /I \ duk 



r-2 + fs + 3(5-e 5 + 



^) 



Uk = 0, (3.28) 



which can be written in the form 



y 



2(fuk 
dy"^ 



+ il-2p)y^+[eV+p'-'^'] Uk=0 



(3.29) 



with 



P = ^(1 + s), 
i={l-e-sr\ 

z.= 3+ s- 5 + 3efl + ^ 



(3.30) 
(3.31) 
(3.32) 
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Figure 1. For comparison between the numerical and analytical predictions, the value of was 
estimated using (3.38) for 13056 parameter sets with V-r, in the observed range. The histogram shows 
the percentage error between the estimated result and that obtained from numerical integration of 
the perturbation equations. It can be seen that the vast majority of the ris-values estimated using 
(3.38) differ less than 1% from the calculated value. 



The solution of eqn. (3.29) is Uk = z^(i^{iz) where is a Bessel function of order i' [32]. 
The general solution is 



Uk{y) = y2+§ \cM^\iy)+C2HP{£y) 



(3.33) 



using Hankel functions of the first and second kind. Comparing with the short wavelength 
solution Uk{y » 1) ~ e"^'^"'^ I \j2csk means C2 = and c\ = y^vr/lAJcJ so (up to a phase 
factor) 

' ' ' y 



My) = iwA- 



y 



2V c,/c V 1-e-s " \1 



e — s 



I — 3 

In the long wavelength limit i/^^^(y < 1) ~ J l^^'"'"''^'^''''^ rmky''' ' ^° 



r(^) 



Uk{y)\^2''-2-^(l-e-sr~2 



1 



r(i) 



Thus, from (3.24) we find 



1/2 yj2cskHy 



n 



23/2 



(3.34) 



(3.35) 



(3.36) 



Defining V(i^) = 2'^ ^(1 — e — sY 2r(t')/r(|) and using (3.35) this can be written 



V. 



1/2 (V{u)\ H 



n 



vr 



/Coe 



so the spectral index is 



n, = 1 + (3 - 2v), 

= 1 - 2s + 2(5 - 6e f 1 + ^ 



(3.37) 



(3.38) 
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This estimate compares extremely well (significantly less than 1% error in most cases) with 
the value obtained by solving the perturbation equations numerically described in the fol- 
lowing section. The comparison is illustrated in Fig. 1. Note that, as (3/r] is not necessarily 
less than one, the terms in ef^/rj are of the same order as the other slow-roll parameters and 
should not be neglected in (3.38). In fact, numerically, we find that dropping these terms 
shifts the value of the spectral index by a large factor of roughly 5-10%. Comparing (3.38) 
with the corresponding equation in single field slow-roll Ug = 1 + 25 — Qe [33] and the single 
field DBI case^ Ug = 1 — 2s + 25 — [7], it can be seen that the standard expression is 
recovered when the coupling is small. 

In Fig. 1, (3.38) was evaluated at sound horizon crossing kcs = aH i.e y = I. As the 
definitions of e and 5 explicitly include factors of ^ = e^^ , it is useful to check to what extent 
the power spectrum and the spectral index depends on the evolution of the fields. The sound 
horizon crossing formalism can be expressed by the condition 



d f H 3 

— — — y"^' 

dy Vvcle 



0, 



(3.39) 



1/2 

[30] which ensures that the that the expression for V-j^ is independent of y and can safely 
be evaluated at the sound horizon crossing y = 1. Following [30], we can generalise this 
condition. Observing that 

d_ _ 1 d 

dy 



yH{l 



df 



(3.40) 



it follows that 



d 



^ dy 



d\og{y) 



log 



H 3 



H 

\fcle 



3 



H 



V + 



/c^e H{1 — e — s 
e + y + e/{2eH) 
1-e-s 

-5 + 3e(l + ^ 



1 



d 



dt 
H 



H 



/c,,e 



3 

:y2^ 



3 

y2" 



(3.41) 



Inserting the value of v in (3.32), obtained by neglecting the two-field correction to the Uk 
mass term, it can be seen that (to first order in the slow-roll parameters) the expression for 
the power spectrum is independent of the y-value at which it is evaluated. This justifies our 
assumption that we can evaluate the power spectrum at sound horizon crossing. 



4 The parameter space of the model 

The scalar-tensor DBI inflation model has thus far been investigated in detail in scenarios in 
which the 99-field is in the minimum of its effective potential, in order to specify uniquely the 
initial conditions and to neglect short-term effects due to the field oscillating as it reaches 
the minimum^. Even in this restricted case, the inflationary dynamics are determined by six 
parameters: the coupling the t'Hooft coupling A and the cutoff scale /x in the DBI warp 

^ This result is quoted [cf. eqn. (40) in [7]] in the form Us — — 2e — fj — s where fj = e/{He). However, one 
can see from (3.11) that in the limit /3 — 0, 77 = s — 25 + 4e, which gives the stated result. 

^ In this case the initial value of the x-field is of less importance as, for a given set of parameters, the 
system proceeds along a single trajectory in field space. The anisotropies in the CMB radiation are affected 
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Figure 2. Spectral index plots for the 'better' set (with power spectrum amplitude within obser- 
vational limits) for the exponential potential, colour-coded according to the value of 11^. Upper- left 
panel: the j3 — logj^g m plane; upper-right panel: the logj^g Uq — rj plane; lower-left panel: logj^g \ — n 
plane; lower-right panel: logj^g m — logj^g Uq plane. 

factor, the mass term in the quadratic potential of the DBI field m and the two constants 
in the (/j-field potential Uq and r]. It is possible to understand the importance of these 
different parameters by keeping five fixed and varying one (as in figs. 5, 7 and 9 in [27]). 
However, this only gives a limited view of the effect of changing a particular parameter on 
the observations. To improve upon this, one can perform multiple runs using random sets of 
parameters from within a given range to compare the resulting values of the power spectrum 
amplitude and spectral index with observations. In the following, the results of such an 
analysis are analysed for the exponential potential. The quadratic potential yielded similar 
results: these are discussed in Appendix B. 



by the range of scales (around 4 orders of magnitude or 10 efolds) that cross the curvature horizon roughly 55 
efolds before the end of inflation. Thus, only the behaviour of the system in the last 60-70 efolds of inflation is 
needed to compare with observations and we are free to choose values Xini that allow an inflationary period of 
sufficient duration, bearing in mind string theory constraints on the maximum length of the throat (cf. [5, 7]) 
that mean Xini ^ in Planck units. In all of the examples presented here, Xini = 1-5 has been used. 
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(7 - l)|fc = f.,ff 7max- 1 



Figure 3. Correlations between tire spectral index and the DBI boost factor for the exponential 
potential (Note the log scale on the x-axis.) The left panels show (7 — 1) (calculated at fcpivot — aH) 
against n^, with the green region indicating the observational limits on the spectral index. The right- 
panels show the correlation between (7max — 1) and rij. (Points are colour-coded according to the 
value of Tig, for easy comparison with the parameter-space plots in figs. 2). 



The parameter ranges^ used in the run with the exponential potential were 

/5 G [0.1, 4], A G [lO^O-S, 5 X 10^3], m G [lO'^-^, lO'^'^], 
Uo G [10-^■^ 10-6], fi G [0.01, 0.5], r] G [0.5, 3.5]. (4.1) 

5 million parameter sets were investigated, of which 2,231,657 (44.6%) satisfied the criteria 
for the background evolution^. Most of the points rejected at this stage (60%) yielded too 
few efolds, with the majority of the remainder exhibiting slow-roll in x (38%). With the 
range of parameters given above, the volume of the parameter space was rather large, which 
meant that only a relatively small number of parameter sets (13063) were found that gave 
rise to a power spectrum amplitude within the observed range^. This shall be referred to 

^ These ranges were based on preliminary tests in which it was found that a large range of values of the 
parameters Uo and A (and to some extent, m) gave rise to a power spectrum amplitude within observable 
limits. Only non-zero couplings /3 > were considered so as to ensure the presence of a minimum for the 
</p-field, in which case there is a unique trajectory in field space for a given set of parameters. 

® The two most important criteria are that there should be a sufficient amount of inflation to solve the 
horizon and flatness problems and that the maximum value of 7 (considered over the entire run) should be 
greater than 1. The latter case, of course, is far from unphysical as it corresponds to a two-field model with 
slowly rolling coupled scalar fields. However, in this section we focus on deviations from the standard scenario 
so exclude these models to focus on the more unusual phenomenology associated with the DBI field. 

This can be understand by considering the analogous process in the single-field case, in which the power 
spectrum amplitude Pamp is proportional to the squared mass of the field. Increasing the volume of the 
parameter space (either by extending the range of mass values tested or increasing the dimension of the 
parameter space) means that the number of cases that do not give rise to the observational value of Pamp is 
vastly increased. 
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Figure 4. The /3 — log^o Uq (left column) and logj^g A — logj^g ™ (right column) planes for the 'better 
set' for the exponential potential. The upper plots are colour-coded according to the number of cfolds 
in the run iVmax and the lower plots according to the value of A at fcpivot = aH. 



as the 'better' set. Similarly, the 'best' set for which both the power spectrum and spectral 
index were within observational limits yielded 3834 points. 

Except for some runs with the coupling f3 very close to zero, no distinct region of the 
parameter space was excluded by discounting the parameter sets that did not satisfy the 
criteria for the background evolution (i.e. gave rise to an extremely short inflationary stage 
or standard slow-roll models). The relatively large range of parameters considered means that 
the calculated value of the power spectrum amplitude varies over several orders of magnitude 
with a strong dependence on the mass scale of the DBI field (as is to be expected since it 
is the perturbations of the DBI field that are the dominant contribution to the curvature 
perturbation). Larger values of the (/J- field potential, corresponding to small r] and large Uq, 
increase the Hubble parameter and also give rise to larger values of Pamp- 

Focusing on the 'better set' of runs with Pamp within observational limits, one can see 
in fig. 2 that the parameter space is significantly constrained compared to the 'good set', 
although the range of values of the spectral index is quite large. There appear to be 
two overlapping regions in the plots in fig. 2. The light blue points with ~ 0.97 make 
up the bulk of the 'best set' (shown in more detail in fig. 5) and have 7 ~ 1 as the pivot 
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scale crosses the horizon, so the spectral index is similar to that in slow-roll models as DBI 
corrections are minimal. To illustrate more clearly the qualitative difference in the properties 
of the background runs that give rise to the distinct regions in fig. 2, correlations between 
the boost factor and Ug are shown in fig. 3. As 7 is a dynamical quantity that increases to 
a maximum as inflation proceeds, two views are shown: in the left panel, the value taken 
when the pivot mode leaves the horizon (7 — l)|fc ^^^^^^^^ (corresponding to = 55) and the 
maximum value obtained throughout the background run (7max — 1) (right panel). As the 
value of Hg shown on the vertical axis is a constant property for a given run, the points can 
be thought of moving horizontally between the left panel and the right as inflation precedes. 
It can be seen that the runs which do exhibit significant DBI characteristics (as indicated 
by 7max ^ 100 in the right-panel plot) are precisely those with larger values of the boost 
factor at horizon crossing that are excluded by virtue of a blue-tilted spectral index. This 
is, however, not to say that there are no viable scenarios that exhibit DBI characteristics 
as some of the darker blue points — which, as can be seen by comparing figs. 2 and 3, are 
associated with smaller values of A and ^, and a DBI mass scale m ~ 0(10"^) — correspond 
to runs in which the boost factor increases from a very tiny value to 7max ~ 0(10). In 
both cases, deviations from the quasi-slow-roll solutions with 7—1^1 throughout the 
background evolution give rise to a corresponding shift in the resulting spectral index of 
the perturbations, leading to either a blue-tilted (for the points corresponding to solutions 
bearing close resemblance to standard DBI models) or a more red-tilted spectrum (for the 
points in which DBI characteristics are suppressed by the presence of the second field). 

The parameters A and m strongly affect the other background quantities: the duration 
of inflation A'max and the coupling A (measured at /cpivot = aH). The models with A > 10 
(and Pamp in the observed range) have small m to compensate (cf. fig. 4). In a similar way 
to the case shown in fig. 9 of [27] these models can exhibit slow-roll inflation in x so that 
A^max is very large. As might be expected, models with large A have /3 > 1; the left column 
of fig. 4 shows a noticeable correlation between f3 and logj^g ^0 for these models. 

Looking at the 'best' set for this model, with both rig and Pamp within observational 
limits (fig. 5) it can be seen that most of the allowed models have very similar properties, 
with the boost factor 7 at horizon crossing deviating only slightly from 1 and rising to values 
less than 0(10). The allowed range of the DBI mass m is relatively small (one order of 
magnitude) compared to Uq and a large range of /3 values are compatible with observations. 
The plots showing the (3 — log^^g ^0 and log^g — rj display clear structure, favouring larger 
values for the potential of the (/j-field (when /3 is not large) than in typical slow-roll scenarios, 
further emphasising the non-dynamical role of the field. The bottom-right plot, showing the 
log]^o — plane is also tightly constrained, with values that lead to large f{x) excluded. As 
in the standard DBI scenario, a larger value of the warp factor is related to a large number 
of efolds of inflation. 

5 Conclusions 

In this article we have investigated in detail the consequences of the scalar-tensor model of 
DBI inflation introduced in [27], paying particular attention to the effect of the different 
parameters on the value of the spectral index. When the coupling between the two fields 
is non-zero, the (^-field (with a canonical kinetic term in its Lagrangian) is forced into the 
minimum of its effective potential. Although it does not play a dynamical role, the field 
indirectly affects the DBI field x and the perturbations via the conformal coupling factor 
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Figure 5. Views of tiie parameter space for tlie 'best set' (for wiiicii tiie calculated values of Ug and 
-Pamp are within observational limits) for the exponential potential. Top row: the /3 — logj^Q m (left) 
and logiQ A — logj^g m (right) planes, colour-coded according to the value of log]^o(7~ 1) at Apivot = aH. 
Middle row: the /3 — logj^Q Uq (left) and logj^Q m — logj^Q Uq (right) planes, colour-coded according to 
the value of log]^Q(7niax ~ !)• The bottom row shows the logj^g Uq — r/ plane, colour-coded according 
to the value of A at fcpivot ~ o,H, and the logj^Q X — fi plane, colour-coded according to the value of 



A = exp{f3ip). With the (^-field in its minimum and working in the 'slow-roll limit' in which 
the equation of motion for the DBI field is dominated by the potential V{x) (with effective 
mass uiA) we have shown that, with appropriate slow-roll parameters, the perturbation 
equations can be put in a form similar to that describing perturbations in single-field k- 
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inflation, with an additional small correction due to the two-field nature of the full system. 
In the limit that this can be neglected, the resulting estimate for the spectral index is 

= 1 - 2s + 25 - 6e (l + f ) , 

which is superficially similar to the standard expressions for the single field DBI scenario. 
However, the the e and 6 parameters are proportional to A~'^ (1 + 4/3/r/)~^ and A~'^ (1 + 4/3/r/) 
respectively, so tend to be smaller than their counterparts in the uncoupled case (with all 
other parameters unchanged). The values of the slow-roll parameters, are dependent on the 
background evolution of the fields (which in turn are determined by quantities such as the 
mass of the DBI field) and the actual value of the can exhibit either a blue- or red-tilt. 

Understanding the effect of the free parameters of the model on the observable quantities 
is essential: the numerical work described in Sec. 4 represents a first step toward this end. It 
is shown that the parameter space can be severely constrained by the requirement that the 
mechanism must be capable of driving a sufficient amount of inflation and by considering the 
values of the power spectrum amplitude and the spectral index. The latter measurement in 
particular rules out a large class of models — those with small values of Cs at horizon crossing, 
which give rise to a blue-tilted spectrum — so that only scenarios in which DBI characteristics 
are suppressed at the beginning of the observable period of inflation are allowed. In general, 
the excluded models with (7 — 1) ^ C(l) at horizon crossing exhibit significant DBI effects 
(i.e. large values of the boost factor)throughout the inflationary period. The remaining 
realisations of the model have (7 — 1) <C 1 at horizon crossing: some of which differ only 
slightly from standard slow-roll solutions. This is in general agreement with the results of 
other numerical studies investigating the standard DBI model [5] where it was found that 
only models exhibiting small DBI effects were viable. However, a number of viable cases — 
associated with a more red-tilted spectral index — were found for which the boost factor 
is initially suppressed by the effect of the coupling between the fields, but increases later to 
moderate values ~ 0(10). Further work is needed to fully understand the viability of the 
model; a natural extension is to consider other observables: in particular, the tensor-to-scalar 
ratio r and the level of non-Gaussianity (proportional to 7^ in the single-field case) which 
so tightly constrains the standard DBI scenario. It would also be of interest to calculate the 
running of the spectral index, as including this gives different bounds on the spectral index 
[34]. 

In conclusion, the scalar-tensor DBI inflation scenario provides a simple mechanism to 
reduce the large values of the boost factor associated with single field models with DBI action, 
whilst still being able to drive 60 efolds of inflation. We have indicated the regions of the 
parameter space of the model capable of giving rise to a power spectrum with amplitude and 
spectral index within the observed bounds, and moreover shown analytically that the value 
of the latter quantity arises when the coupling forces one field into its minimum, leaving the 
DBI field to slow-roh. 

Acknowledgments 

D.F.M. thanks the Research Council of Norway FRINAT grant 197251/V30. D.F.M. is 
also partially supported by project PTDC/FIS/111725/2009 and CERN/FP/116398/2010. 
J.M.W. was supported by Yggdrasil Grant 195757/Vll from the research council of Norway, 
and would like to thank the members of the Institute of Theoretical Astrophysics at the 



- 16 - 



University of Norway for their kind hospitality while this work was initiated. We would also 
like to thank the anonymous referee for helpful and constructive comments. 



A Approximating C^^ to first order in slow-roll parameters 

In this appendix, we derive an first-order expression for the mass term Cy.y. [given in full 
by (3.15)] that appears in the perturbed equation of motion (3.14). Before rewriting each 
contribution to this equation in terms of the slow roll parameters e, 5 and s, it is necessary 
to obtain an expression involving To do this, one can differentiate (1 — c^) = A~'^fx^ to 
get a relation between x and This yields 



1 fx. 



^ =/3vi-|^X+(l 
X 



> r. ' 



X 
Hie 



^^_ifxX 



+ 1 



(A.l) 
(A.2) 



Then differentiate —2H ~ A'^jx^ to obtain, 



2€H^ + ieHH ~ {2/3ip - c./c,) {2€H^) + 2A^^XX, 
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Substituting (3.11) and (A.2) gives 

_ iIk^ r 
' f H- 



(l-c^) \-5+(l + 



213 



e. 



(A.4) 



We can now write each contribution to C^^ in (3.15) in terms of the slow roll parameters e, 
6 and s, working to 0{e) and dropping second order terms. Starting at the beginning 



The first term is 
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using A'^jx^ ~ 2eif^. 

• The second term is 
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The third term is 

'ifxX 



1 fx (1 „2\t a2. 



3H^ 



(A.5) 



where (3.7) was used in the third step and (A. 4) in the last. 
The fourth term is 



2, +{l + Cs)f 



-1 I fx 

f 



2(1 + c. 



(1 + 2Cs) - Cs 



The final part of this term can be expressed in terms of time derivatives of the slow-roll 
parameters and other second-order terms 



1 (f. 
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(A.6) 



using (A. 4) and (A. 2). We have checked numerically to confirm that the contribution 
of this term to C^^ is negligible. 



The fifth term is 



The sixth term is 



The seventh term is 



UW(l + c^)-3i/^(l + c^)e. 



(A.7) 
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using (3.10). 
• The eighth term is 
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Figure 6. The distribution of (7 — 1) values for the offset quadratic potential. Quantities plotted 
are the same as in fig. 3. 



The final term is 



XX 



XX 

V 



(A.9) 



Combining eqns. (A. 5), (A. 7), (A. 8) and (A.9) then gives eqn. (3.16) 



Cj^^ 2^ 3HU S - s - 2e 



B Comparing the exponential and quadratic potentials 

As mentioned in the body of the article, the numerical analysis in Sec. 4 was repeated with 
an offset quadratic potential for field (as described in [27]) given by 



Uiip) = Uo{v-r^f. 



(B.l) 



Like the exponential potential, this is dependent on two free parameters Uq and rj: in this 
case the latter gives the value of the true minimum is (p = t]. The parameter ranges were 



/3 G [0.1, 4], A G [10lO■^ 5 x 10^^], m G [10-^■^ lO"^-^] 



Uo G [10~^\ 10 



fi G [0.01, 0.5], r] G [1, 6]. 



(B.2) 



Again, 5 million parameter sets were investigated, of which 2,748,236 (55.0%) satisfied the 
criteria for the background evolution. As with the exponential potential, of the points rejected 
at this stage, most (56%) yielded too few efolds, with the majority of the remainder exhibiting 
slow-roll in x (41%). 11,999 of the acceptable parameter sets had a power spectrum amplitude 
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Figure 7. Two views of tlic parameter space for tlie 'best set' (for wlriclr tlic calculated values of Us 
and Pamp are within observational limits) for the quadratic potential, colour-coded according to the 
value of A at fcpivot = ciH. Left panel: logj^Q A — /i; right panel: logj^Q m — log^Q Uq. 

within the observed range, and the 'best' set (for which both the power spectrum and spectral 
index were within observational limits) yielded 2548 points. 

The parameter space for the quadratic potential is very similar to that of the exponential 
potential in many respects, with the exception of those involving the parameter r], which has 
a different interpretation in each case. For example, for the exponential potential smaller 
values of r] give rise to larger values of the power spectrum amplitude, due to the increase in 
size of the potential U{ip). For the quadratic potential, for the same reason, it is larger values 
of 7] that are correlated with Pamp- Focusing on the parameter space for the 'better' set for 
this potential consisting of parameters that give rise to -Pamp within observational limits, it 
is again possible to identify two groups of points in the parameter space. The interpretation 
is the same as in the case of the exponential potential: runs that deviate significantly from 
quasi-slow-roll behaviour (i.e. (7 — 1) ^ 1) as the pivot mode leaves the horizon give rise to 
a blue-tilted spectrum with > 1 (blue points) and runs in which DBI characteristics are 
suppressed at = 55, but become significant later during inflation, exhibit a strong red-tilt 
(dark red points). The similarity between the two potentials can be seen by comparing figs. 
3 and 6, which show the correlation between the spectral index and the boost factor. 

Looking at the 'best' set, with both the power spectrum amplitude and spectral in- 
dex satisfying observational constraints, we see more similarities between the results for the 
quadratic and exponential potentials. The distribution of points in the views of the param- 
eter space showing the mass parameters and the warp factor parameters shown in fig. 7 are 
almost identical to those for the exponential potential (cf. middle-right and bottom-right 
plots in fig. 5). Again, there is a strong correlation between small values of the DBI mass 
and large values of the coupling A (and also large values of A^max)- 

The general distribution of the points in the parameter space is largely the same for 
both potentials although the (3 dependence is a little different. For the exponential potential, 
the set of allowed parameters is cut off abruptly in the /3 — logj^Q Uq plane (see in the left- 
hand plots in fig. 4) when (3 and Uq are large as the Hubble scale (and by extension, Pamp) 
becomes too large. In the equivalent plot for the quadratic potential (not shown) H depends 
also on the displacement of the field from its true minimum, and so the effect of Uq and /3 is 
modulated by other parameters and the cutoff is not so clear. 
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Figure 8. The (3 — rj planes for the quadratic potential (left) and the exponential potential (right), 
using the 'best sets'. Both plots are colour- coded according to the value of iVmax (cf. analytical 
constraints shown in fig. 9). 




P P 

Figure 9. For comparison with fig. 8, the shaded regions show the regions of the rj — /3 planes in 
which the condition pdbi > P:p is satisfied for the quadratic (left panel) and exponential potentials 
(right panel) using the same parameter ranges. The former condition is dependent on the ratio V/Uq 
(see text) so we have used V/Uq = 10"^, 10"*, lO"-^, 10~^, with the lightest shade corresponding to 
10"^ 

We have remarked that the differences between the observationally viable regions of 
the parameter spaces for the exponential and quadratic potentials are manifest when the 
parameters /3 and r] are plotted. It is interesting to compare the plots in fig. 7 with those in 
fig. 9, representing the constraints derived analytically (in [27]) under the assumption that 
the fields are dominated by their potential terms. As shown in [27] the condition poBi > Pip 
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is given by W{x ~ /Uq) < 2 in terms of Lambert's W function for the quadratic 

potential and rj > 4f3 for the exponential potential. It can be seen that the areas in fig. 8 that 
correspond to the shaded regions in fig. 9 — in which the DBI energy density is dominant — 
are almost completely devoid of points. Also, many of the points are clustered in the vicinity 
of the pdbi = P-fi boundary. Thus, in the majority of runs with parameters that give rise to 
a power spectrum with a spectral index within observable limits, the energy density of the 
DBI field is a subdominant but non-negligible fraction of the total energy density. 
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